On the interaction of the electromagnetic radiation with the breaking plasma waves 
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An electromagnetic wave (EMW) interacting with the moving singularity of the charged particle 
flux undergoes the reflection and absorption as well as frequency change due to Doppler effect and 
nonlinearity. The singularity corresponding to a caustic in plasma flow with inhomogeneous velocity 
can arise during the breaking of the finite amplitude Langmuir waves due to nonlinear effects. A 
systematic analysis of the wave-breaking regimes and caustics formation is presented and the EMW 
reflection coefficients are calculated. 
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I. INTRODUCTION 

A strong interest has persisted in studying the elec- 
tromagnetic wave (EMW) interaction with the relativis- 
tic electron structures represented by electron beams [1] , 
Langmuir waves [21 |3] H] [5j [6] ionization fronts [7] , rela- 
tivistic solitons and vortices [5] in an undcrdcnsc plasma, 
and by the oscillating electron layers at a solid target sur- 
face [amain]. 

The laser produced irradiance approaches 10 22 W/cm 2 
12J, which leads to ultrarelativistic plasma dynamics. By 
further increasing the irradiance we shall see novel physi- 
cal regimes such as the radiation friction force dominated 
EMW-matter interaction [13]. At irradiances around 
10 28 W/cm 2 [15], the focused light becomes so strong 
that the nonlinear effects come into play, as predicted by 
quantum electrodynamics, including the vacuum polar- 
ization and the electron-positron pair creation from vac- 
uum (see review articles [16 and the liteterature quoted 
therein) . The studies of the EMW interaction with the 
relativistic structures pave a way towards achieving such 
the irradiance. 

Here we pay the main attention to the "flying mirror" 
concept [4] (see its further development and discussion in 
Refs. [IU HE1 EZ] ) • The electromagnetic radiation source 
suggested in this concept has unique advantages. It is 
robust since it is based on the fundamental process of 
wave breaking. It enables the control of the output pulse 
parameters such as frequency, duration, focusing, spec- 
trum, etc., required by a wide range of applications. It 
allows so high up-shifting of the laser radiation frequency, 
followed by focusing to a spot whith size determined by 
the shortened wavelength, that the quantum electrody- 
namics critical field (Schwinger limit) can be achieved 
with present-day laser systems [4]. 

In the "flying mirror" concept, extremely high den- 
sity thin electron shells in the strongly nonlinear plasma 
wave, generated in the wake of an ultra-short laser pulse, 
act as semi-transparent mirrors flying with a velocity 
close to the speed of light, Fig. [l] The mirrors reflect a 
counter-propagating pulse, whose frequency is upshiftcd, 
duration is shortened and power is increased. The ob- 




FIG. 1: (color online) The "flying mirror" concept. The pulse 
incident on the moving density cusp is partially reflected. The 
reflected pulse frequency is upshifted. 



taincd pulse power is proportional to the Lorentz factor 
jph = (1 — 0tb)~ X corresponding to the velocity of the 
mirrors v p h — f3 P h,c while the frequency and the pulse 
compression are proportional to the square of the Lorentz 
factor. 

A nonlinear wake wave with substantially large am- 
plitude, at which the Lorentz factor of electrons, j e , is 
greater than that of the wake wave, "f e > "f p h, breaks 
and forms cusps in the electron density distribution |19j . 
It is the steepness of the cusps that affords the efficient 
reflection of a portion of the counter-propagating laser 
pulse. If the wake wave is far below the wave-breaking 
threshold, the reflection is exponentially small. The fre- 
quency of the reflected light is upshifted due to the double 
Doppler effect, as predicted by Einstein [20 . In his pa- 
per, the problem of the reflection of light from a mirror 
moving with constant velocity, vm — C0M> close to the 
speed of light , is solved as an example of the use of rel- 



ativistic transformations. The incident and the reflected 
wave frequencies, ujq and u> r , are related to each other as 



u> r = ujq- 



l + l3 2 M + 2p M cos6 



1-/3 2 



(1) 



A I 



where 9$ is the wave incidence angle. The reflection an- 
gle, 9 r , is defined by the expression 



cos6L 



(l + /?^)cosflo + 2/? M 
1 + 2(3 M cos 6 + [3 M ' 



The wave amplitude transforms according to 



E ,. 



R^En 
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(3) 



where R is the reflection coefficient. The role of the 
moving mirror is played by a steep cusps of the elec- 
tron density in the wake wave, therefore we write (3m — 

Pph = v ph /c and -f M = l P h = Va/ 1 ~ P 2 p h- Tne wave 
propagates in a narrow angle A8 jw 1/jm- As a result 
of the interaction with the moving mirror, the reflected 
electromagnetic pulse is compressed and its frequency is 
upshifted by a factor which in an ultrarelativistic limit, 
Iph 3> 1, is approximately equal to 47^ f . It is also con- 
ceivable that the interaction of the counter-propagating 
pulse with the thin relativistic electron shell leads to a 
higher order harmonic generation, according to a mecha- 
nism considered in details in Refs. [5J[TB]. The harmon- 
ics allow additional compression with the factor ~ An^f,j, 
where n is a harmonic number. 

It is important that the relativistic dependence of the 
Langmuir frequency on the wave amplitude results in the 
formation of wake waves with curved fronts that have a 
form close to a paraboloid [TS] . The paraboloidal mirror 
focuses the reflected pulse, which results in the intensi- 
fication of the radiation. In the reference frame moving 
with the mirror velocity the reflected light has the wave- 
length equal to rs A /2jm- It can be focused into the 
spot with the transverse size Xo/2^/m with the intensity 
in the focus given by 



R/o7m 



(4) 



Here Iq is the incident wave intensity at the pulse waist 
D. 

The flying mirror formed by an electron density mod- 
ulations in the breaking wake wave can also transform 
the quasistationary electromagnetic field into a high- 
frequency short electromagnetic pulse [5], such as a low- 
frequency electromagnetic field of a relativistic electro- 
magnetic sub-cycle soliton, quasistatic magnetic field of 
an electron vortex, and the quasistationary longitudinal 
electric field of a wake wave. 

The "flying mirror" concept was demonstrated in the 
proof-of-principle experiments [5) |B], where the narrow- 
band XUV generation was detected. This opens a way for 



developing a compact tunable coherent monochromatic 
X-ray source with the parameters required for various 
applications in biology, medicine, spectroscopy and ma- 
terial sciences [6 . 

Reasoning from the principal importance of the non- 
linear wake wave dynamics in the "flying mirror" con- 
cept, in the present paper we consider a variety of the 
wake wave breaking regimes. We find typical singulari- 
ties of the electron density distribution in the breaking 
wave and calculate their reflection coefficients, R. This 
provides a key for optimising an operation of the hard 
electromagnetic radiation source based on the ' 'flying 
mirror" concept. We describe the EMW interaction with 
the multiple cusp structure and show the enhancement of 
the efficiency of the light reflection in this configuration. 



II. GENERIC PROPERTIES OF THE WAKE 
WAVE BREAKING 

A. Gradient catastrophe 

The wave breaking in collisionless plasmas provides an 
example of typical behaviour of the waves in nonlinear 
systems. The most fundamental properties of nonlinear 
waves can be represented by Riemann waves in gas dy- 
namics (e.g. see Refs. [3TJ[22]), described by the equation 



dtU + w{u)d x u = 0, 



(5) 



where u{x, t) is the gas velocity and w(u) = u + c s (u). 
The sound speed c s (u) is a given function of u, in a par- 
ticular case of a cold gas of nonintcracting particles it is 
equal to zero. The solution to Eq. (I5J) can be written in 
an implicit form 



u = Uq(x — w(u)t) 



(6) 



with Uo(x) being an atrbitrary function deternmined by 
the initial condition: u(x)\ t=0 = Uq(x). In a finite time 
the Riemann wave breaks at the location where the gra- 
dient of the function u[x 1 1) becomes infinite. Taking the 
derivative of the expression (l6l with respect to x, we ob- 
tain 



d T u = 



U' {x - w(u))t) 
l + U^(x-w(u)t)w'{u)t' 



(7) 



Here U' and w' stand for the derivatives of the functions 
Uq and w with respect to their arguments. For non- 
trivial dependences of Uq(x) and c s (u), we find that the 
denominator in Eq. (17| vanishes at some point xt, r , i- c. 
the gradient u(x,t) tends to infinity, at time tf, r equal to 
tbr = —l/Ub(xbr ~ w(u)tb r )w' (u) while the gas velocity 
U(, r = Uo(xb r — u> / (w& T .)t& r ) remains constant. This phe- 
nomenon is known as the "gradient catastrophe" or the 
"wave breaking" . 

Another approach for description of the wave breaking 
uses a perturbation theory [22j to find the solution to the 



equation (pH). We write 



u = u (0) + eu^ + e 2 u^ + 



(8) 



with £ <C 1. We assume that in zeroth order the wave 



amplitude is homogeneous with the velocity, 



,(0) 



stant in space and time, 
amplitude, we have 



To the first order in the wave 



d t u {1) + w(u^)d x u (1) = 0. 



(9) 



This is the simplest wave equation which describes the 
wave with the frequency and wavenumber related to each 
other via the dispersion equation u) = kw(u^ '). Thus 
we obtain the wave propagating in nondispersive media 
where both the phase velocity, ui/k, and the group ve- 
locity, dui/dk, are equal to w(u^). The solution to Eq. 
(fOJ) is an arbitrary function of the variable x — wot, where 
Wq = w(u^). We choose it in the form 



r '(X) — u m sink(x — Wot). 



(10) 



To the second order in the wave amplitude, we obtain 

d t u^ - w{u^)d x u^ = -«««/(«W)5 x «W = 

bk 
= -—sin(2k(x-w t)) (11) 

with b = u^w' (u^) . The solution to this equation, 



,( 2 ) 



bkt 
(x,t) = -—-sin(2k(x-w Q t)), (12) 



describes the second harmonic with the resonant growth 
of the amplitude in time. 

To the third order in the wave amplitude, we can find 
that, in general, the amplitude of the the third harmonic 
grows with time as t 2 , and so on. In the media without 
dispersion high harmonics are always in resonance with 
the first harmonic. The resonance between harmonics ap- 
pears because of the fact that the velocity of propagation 
is the same for all harmonics; it does not depend on the 
wave number. This leads to the increase of the velocity 
gradient (the wave steepcnning) linearly with time and 
to the break of even weak but finite amplitude wave. 

The situation with the nonlinear Langmuir wave break- 
ing is different. For simplicity we assume that the wave 
amplitude is nonrelativistic and the ions are immobile. 
We cast the equations of the electron fluid motion and 
the electric field induced in the plasma in the form 



dtv + vd x v 



m, 



-E, 



dtE + vd x E — 4TTeno(ex)v. 



(13) 



(14) 



Here no(ex) with s <C 1 is the ion density which is as- 
sumed to be weakly inhomogeneous. Although with the 
use of the Lagrange coordinates a solution to Eqs. ( 13|14 ) 



can be reduced to quadratures, we use a perturbation ap- 
proach in order to analyse whether or not the high har- 
monics are in resonance with the first harmonic in the 
case of nonlinear Langmuir waves. We expand v(x,t), 
E(x,t), and no{sx) into series: 



v = ev^ ' 

E = eE (1 



+ £ V 2 > + .. 



; 2£(2) 



(15) 



In first order from Eqs ( 13|14 l we obtain the equations 

(16) 



dtV lX> + _^£(D = o, 



d t E {l) -47ren (0)w (1) =0 
with the solution 

v — "m sm(u! pe t — kx), 



E (D 



^^e^m^pe 



cos(cu pe t — kx), 



(17) 



(18) 



(19) 



where lo p& = ^/47rno(0)e 2 /m e is the Langmuir frequency. 
The second order in the wave amplitude yields 

d t v {2) + — £P) = ^ sin (2(oj pe t - kx)) , (20) 
m e 2 



d t E {2) -47re7i (0)w (2) = 

= me ^ efe [1 - cos (2(w pe t - kx))] 



+4iren' Q (0)xv m sin(a; pe i — kx). 
We obtain for v^ and E^ 



(21) 



,(2) 



£[i + -2 ( „ P « t -M]+"i^ 



pc 



x [2to pe t cos (uj pe t — kx) — sm(tu pe t — kx)] , (22) 



E& = 



^ feme -sin(2(^-M)+ Wo(0)OT " 



2e 
x [2ujp e tsin(oj pe t — kx) — cos(uJ pe t — kx)) . (23) 

As we see, in the homogeneous plasma, n — const, 
there is no a resonance between the modes in the Lang- 
muir wave, in contrast to the Ricmann waves consid- 
ered above. This is due to the difference between the 
group and the phase velocity of the Langmuir wave: its 
group velocity in the cold plasma is equal to zero while 
the phase velocity is finite and is given by the relation 
v p h = uj pe /k. The Langmuir wave break occurs if the 
wave is excited with such a strong amplitude, vq, that 
^o > v ph, out of the applicability of the approximation 
used in Eqs. ( |13|14[ ). As seen in Eqs. ( |22|23[ ), the break- 
ing occurs also in inhomogeous plasma, where due to the 
phase mixing effect the wavenumber, k, grows with time 



B. Structure of the Breaking Relativistic Wake 
Wave 



In the particular application for the light intensifica- 
tion [1], the characteristic features of the electron den- 
sity modulations play the key role in the calculating the 
EMW reflection coefficient. 



1. Wake wave breaking in the relativistic electron beam 

We start a discussion of the nonlinear wake wave from 
the consideration of nonlinear perturbations produced by 
the laser pulse in the relativistic electron beam. We ne- 
glect the space charge effect assuming that the laser pulse 
— beam interaction can be described in the test particle 
approximation. The equation for the electron momentum 
reads 



d t p + vd x p = 



2l 



-d~ 



(24) 



Here the longitudinal component of the electron velocity, 
v, along the laser pulse propagation direction, is related 
to the longitudinal component of the electron momentum 

by v = p/^fm e ; 7 = y 1 + a 2 + (p/m e c) is the Lorentz- 
factor and a = eE$/m e ijjQC is the normalized electric field 
of the EMW. We assume that the EMW is circularly 
polarized and depends on the coordinate and time as 
a(x — Vphi). The phase velocity of the wake, v p h, equals 
the laser pulse group velocity, PSJUIIISH]- We note that 
the EMW group velocity, in general, is below the speed 
of light in vacuum if the laser pulse propagates inside a 
waveguide, in a plasma, or/and in the focus region (e.g. 
see Ref. [3D])- Introducing a new variable 



X = x - v ph t, 



we find the integral of Eq. ( 24 1 



7 



PphP 
m e c 



= h 



(25) 



(26) 



This integral plays a role of the Hamiltonian in terms 
of variables X and p, Fig. [2] with h being a constant 
and /? p h — v h /c. If the laser pulse is of a finite duration 
and the electron momentum is pq before the interaction 
with the laser pulse, then the constant h is given by the 

expression h = yl + (po/m e c) — (3 P hPo I 'rn e c. When the 
electron is interacting with the laser pulse its velocity is 
given by the expression 



h(3 ph - Jh 2 -(l-^ h )(l + a 2 (X}) 

v = c 

h-p phy !h 2 -(l-f3l h )(l + ai(X))' 

The solution of the continuity equation, 
d t n + d x (nv) = 0, 



(27) 



(28) 




FIG 
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2: Contours of the Hamiltonian Eq. 
v p hCi,p)-plane for the Gaussian laser pulse a(X) 



a m exp(^— a~j. Shaded region represents trajectories of parti- 
cles whose initial velocity is in the interval < v < v p h. If the 

trajectory with h — 1 is the separatrix, then a m = v h 2 — 1, 
while an enclosed trajectory means a m > \/h 2 — 1. 



gives for the electron density n — n (v p h — vq)/(v p \ 1 — v), 
which is equivalent to 



VW^l h-f3 ph Jh 2 ~l-a 2 (X) 

n = n a - = . (29) 

h - P P Wh 2 - 1 yjh* - 1 - ePQL) 

Here vq = cpo/(m e ch + f3 p hPo) and 



\A^ 



Pb 



(30) 



At the point X 



(29 1 vanishes, i.e. 



= Xi, r , where the denominator in Eq. 

the electron velocity becomes equal 
to the phase velocity v p h, the electron density tends to 
infinity. This corresponds to the breaking of the wave 
induced by the laser pulse in the electron beam. The 
threshold corresponding to the lowest amplitude of the 
laser pulse, at which the pulse piles up the particles to- 
wars the infinite density, is determined by the expression 



a 2 (X 6r ) = h 2 



1. 



(31) 



In the phase plane of the Hamiltonian system described 
by Eq. p6| , the breaking corresponds to X-points and 
to vertical tangents of the contours of the Hamiltonian, 
Fig. [2] 



First, we consider the case when the condition (31 1 is 



satisfied precisely at the maximum of the laser pulse am- 
plitude, a max = a(Xfe r ). In the vicinty of the maximum, 
the amplitude can be represented by the expansion 

a(X) = a(X fcr ) - ia"(X hr )AX 2 , (32) 

where AX = X — Xf, r — > 0, thus for the electron density, 
velocity and momentum we obtain 



'ph 



and 



P = Pr, 



IAXI 



vi|AX|, 



Pi|AX| 



respectively. Here 



m = -~ 



nohaCK-br) 



Vl 

(3 P hh/ 



(h - f3 ph a(X br )) y / a(X br )a"(X br ) 
c(l - ^ ) A / a ( X|)r ) a "(X 6r )/ft, 



(33) 



(34) 



(35) 



(36) 



J-V, 



1 - ^ and ft = ^a(X 6r )a"(X 6r )/(l - /?* h ). 

The electron motion leading to a pile-up and their 
density increase is shown in Fig. [31 We note that the 



singularity in the electron density distribution ( 33 1 is 



r+AX 



not integrable. The integral J_^£n(X)dX diverges. 
Formally it takes an infinite time in order to form the 
singularity. Behavior of the electron density in the 
vicinity of the breaking point can be described by using 



the continuity equation ( 28 1 , the motion integral ( 26 ) 



and the expansion ( 32 1 . The electron displacement is 



described by the equation 



AX 



Vph 



-vilAXI 



(37) 



where a dot stands for the total derivative with respect 
to time. For the initial condition AX| t= o = AXo the 
solution is 



AX = AX exp(-t/r), 



(38) 



where r = h/{\ — [3 2 h )c^ a(X br )a" '(X br ) is the sin- 
gularity pile-up time. The electron density n — 
no |AXo/AX| at the breaking point grows exponentially: 
n = n exp(i/T). Its growth saturates due to a finite 
number of the electrons in the beam or/and due to the 
space charge effect. The case when the space charge ef- 
fect is taken into account is discussed in the next sub- 
section. For a finite number of the particles, the density 
distribution asimptotically in time can be approximated 

by 



n(X) = noW(X), 
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FIG. 3: Formation of the d ensi ty singularity in the beam of 
particles described by Eq. ( 24 1 for the Gaussian laser pulse 

a(X) = a m exp(— X 2 ), a m = v h 2 — 1, and (3 ph — 0.9. (a) 
Thin lines represent the particle trajectories in the (x, i)-plane 
corresponding to different initial locations of particles in the 
beam. Thick lines stand for the laser pulse envelope at the 
corresponding time, (b) The density of the particles at the 
last moment of time. The dashed line shows the theoretical 
limit Eq. ((33). 



where S(x) is the Dirac delta function, and l b is the initial 
length of the bunch. 

Apparently, no break occurs for the phase velocity 
equal to the speed of light in vacuum because the elec- 
tron velocity is always less than c. In this limiting case, 
/3 p h = 1] w e have h — 70 —po/m e c. The electron Lorentz 
gamma factor, 70, the momentum, po, and the veloc- 
ity, vo, before the interaction with the laser pulse can 
be written via the constant ft- as 70 = (1 + h 2 )/2h, 
Po = m e c(l — h 2 )/2h and v = c(l — h 2 )/(l + h 2 ). Inside 
the laser pulse the electron acquires 7 = 70 + a 2 (X)/2h, 
p = po + m e c a 2 (X)/2h, and the electron density becomes 
equal to 



n = no 



« 2 (X) 
1 + h 2 



(40) 



We see that for the finite amplitude laser pulse the elec- 
tron density remains to be finite. 

Now we consider the case when the laser pulse am- 
plitude is larg er than the threshold determined by the 
condition (31), a max > h 2 — 1. In this case the breaking 



(39) occurs at the slope of the laser pulse envelope, where the 



function a(X) has the expansion 

a(X) = a(X br ) - a'(X fcr )AX. 



(41) 



The local profile of the electron density in the vicinity of 
the point X = Xb r is described by 



nn - *(X) 



(42) 



Here x( x ) is the unit-step Heviside function (x( x ) = 
for x < and x(x) = 1 for x > ) and 



"1/2 



no/ia(Xf, r ) 



(43) 



(ft - /3 ph a(X hr .))V2a(X 6r )a'(X hr .)' 
At the breaking point the electron momentum becomes 



equal to pi/ 2 = fiphhU 1 — /3p h - The particle bounces at 

this point, starting to move in the same direction as the 
laser pulse. Its velocity increases from v p ^ at X = X& r to 



'hP ph + V~h 2 -l 

l + (3 ph Vh 2 -l, 



(44) 



for X — > oo. Asymptotically, the electron density in the 
reflected beam tends to 



h + PphVh? 

n = n o — 



1 



h- (3 p hVh? 



1 



(45) 



If the particles are at the rest before interacting with 



1 and h 



1 



■ft 



then 



the laser pulse, i.e. ,<, — x cum ,<, — A/ , / _l — M . 

the density in the reflected beam is equal to no(l + 
/3p h )/(l — /3p h ). The electron motion demonstrating a 
bounce of electrons in the ponderomotive potential of 
the laser pulse and the corresponding compression of the 
density, as well as the electron pile-up at the breaking 
point, is shown in Fig. 



According to Eq. ( 31 1 the breaking coordinate depends 
on the initial electron momentum. For the electron beam 
with some energy spread it results in a non-trivial elec- 
tron density distribution in the vicinity of the singularity. 
We can easily obtain 



+ OC 



n(X) = ni/2 



x(X-X 6r (po))/(j?o)rfp 

(x-x br ( Po )y/2 



(46) 



with f(po) being the distribution function normalized on 
unity: J_ f(p)dp = 1. We assume that the distribution 
function is given by the "water bag" model, i.e. f(po) — 
x(Po,T - |p |)/2po,T, with p 0tT /m e c < /3 ph . Integration 
in the right hand-side of Eq. (46 1 yields 



n(X) = 



"1/2 
P0,T 



X (X - X 6r (-p ,r))A/X-X 6r (-p ) 
-X(X - X br .(p , T ))v/X-X 6r .(po)l • (47) 
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FIG. 4: The same graphs as in Fig. [3] but for the two times 
greater amplitude of the laser pulse, a m = 2\/h 2 — 1. The 
dashed line in frame (b) shows the main term in the expansion 
of the density, Eq. ((42). 



The above consideration is carried out within the 
framework of the test particle approximation, i.e. ne- 
glecting the space charge effect, applicable when the elec- 
tron density is low enough. We note that the space charge 
effects are negligible in the electron-positron plasmas be- 
cause the electron and positron has equal absolute walues 
of the charge and mass. 

Another important case when one can neglect the 
space charge effects is the limit of extremely high EMW 
intensity in the so-called "radiation pressure dominated 
regime" (for details see Refs. [T31I2UI2S]) briefly consid- 
ered in Appendix A. 



2. The plasma wake wave breaking 

Here we discuss the generic properties of the breaking 
of the Langmuir waves in plasmas. We shall consider the 
case when the wake wave is excited by a high intensity 
laser pulse [26l [27j [28]. This consideration is relevant 
to the breaking of the wake wave excited by the ultra- 
relativistic electron bunch |31j . Within the framework 
of one-dimensional approximation, the relativistic Lang- 
muir wave excitation by the electromagnetic wave can be 
described by the equation for the electron momentum, p, 



d t p + vd x p 



27 



-d x a 



(48) 



In contrast to Eq. ( 24 1 , it incorporates the space charge 



effects via the term with the electric field, E, for which 
we have the equation 



dtE + vd x E — A-KeriQev. 



(49) 



Left in a plasma behind a finite width laser pulse, 
the Langmuir wave has the period T w f = 2n/uj pe = 
y'm e /2ne 2 , in the limit of small amplitude, p m /m e c — > 0. 
In the opposite limit of a large amplitude, p m /m e c ^> 1, 
the wake wave period is T w f — 4y // 2p m /m e c/oj pe . Here 
the maximum wave amplitude, p m , is related to the laser 
pulse amplitude as p m = m e ca 2 /2. Therefore, the wake- 
field wavelength X w f sa cT oc a is proportional to the 
laser pulse amplitude. 

For the laser pulse propagating with constant phase 
velocity v p h, we seek a solution of Eqs. ( 48|49 ) in the 
form of a progressive wave, when the solution depends 
on the variables t and x in the combination ( 25 1 . As a 



result we obtain the ordinary differential equation 



7' 



P 



> r ph 



Upe P 

c 2 7m e v p h - p ' 



(50) 



where the prime denotes a derivative with respect to 
X. For piecewise-constant profiles of a(X), a solution 
of Eq. ( 50 1 can be expressed in terms of elliptic inte- 
grals. Here we analyze the solutions of this equation 
in the vicinity of the singularity. The right hand side 
of Eq. (50 1 becomes singular when the denominator, 
7m e v p h — p, tends to zero, i.e when the electron veloc- 
ity v becomes equal to the pase velocity of the wake 
wave, Vph- In the stationary wake wave, the singular- 
ity is reached at the maximum value of the electron ve- 
locity, v m = p m /jm.m e , at X = X m . This singularity, 
called also the "gradient catastrophe", corresponds to 
the wave breaking. In other words, the Langmuir wave 
breaks when the electric field E w f is above the Akhiczcr- 
Polovin electric field Ea-p = (jn e u} pe cje)^j2(^i p ]- l — 1) 
with 7 p /j = l/i/l — /3ph, which means that the electron 

displacement inside the wave becomes equal to or larger 
than the wavelength of the wake plasma wave (see Ref. 

PI)- 

In order to find the singularity structure we expand 
the electron momentum, p, and the Lorentz factor, 7, in 
the vicinity of their maximum, AX — X — X m — > 0. The 
momentum is represented by 



p = P, n + Sp+0(Sp 2 ) 1 



with 



Pm = m e cP ph 



1 



ff 



(51) 



(52) 



P h 



and \5p/p m \ <C 1. Here a m = a(X m ). Keeping the main 
terms of expansions over 5p in both sides of Eq. ( 50 1 , we 
obtain 



m' 



dp 



(53) 



^)L0 2 vr? e c. Multiplying the left- 



withx 2 = 2/3ph7p h ( 1 - 

and right-hand sides of Eq. ( 53 ) on (Sp 2 )' and integrating 

over X, we obtain 



SpSp' = \Jm\c 2 f 2 — x 2 Sp, 



(54) 



where / is an integration constant. We note that since 
Eq. (53) has singularity due to the initial condition 
<5p| AX _ = 0, in principle, one can choose different val- 
ues for the integration constant in different intervals with 
boundary AX = 0, i. e. / = /~ in the interval AX < 
and / = /+ ^ /- for AX > 0. The behavior of the 
function 5p(X) at the singularity is different depending 
on the behaviour of the product SpSp' at Sp — > 0. 

If the product SpSp' is finite and nonzero for 5p — v 0, 
then the main term in the expansion of the solution of 
Eq. ([54b at AX -> is 



Sp = -m e c^2f\AX\. (55) 

Using this expression we find the electron velocity: 



Vph 



e / 2/|AX| \ 1/2 
t h \ 1 + < ) 



(56) 



The electron density, n(X), is given by the solution to the 
continuity equation ( 28 1 . In the vicinity of the singularity 
it reads 



n ° v P h ~ a 3 ( l + a li \ 



1/2 



Vph 



0, which means / 



(57) 
0, then the 



If at dp -»• 0, SpSp' 
main term in the expansion of the solution of Eq. ( 54 1 
at AX 



reads 



Sp=-[-xAX) 



2/3 



-me/3ph7ph 



■i 1 3uj pe ^/l + al, 
72/3 ph 



2/3 



AX 



For the electron velocity we obtain 



'ph 




2/3 



V2f3 ph 



AX 



(58) 



(59) 



In the vicinity of the singularity the electron density dis- 
tribution is given by 



npVph 
v p h -v 



n "/ph 



'V2(i + C 4) 1/4 A 



2/3 



P h 



3w pe AX 



(60) 



We found three types of the singularity which appears 
in the electron density due to wave breaking. In each case 
the singularity corresponds to the caustic in the plasma 
flow with inhomogencous velocity. General description 
evoking geometric properties of the hydrodynamic flow, 
presented in Appendix B, shows that the electron den- 
sity singularity can be approximated by the dependence 
n(X) tx (X - X m )" Q with 1/2 < a < 1. 



III. ELECTROMAGNETIC WAVE 

REFLECTION AT THE CREST OF THE 

BREAKING WAKE 

As we discussed in Introduction, the "flying mirror" 
concept proposed in Ref. [4] makes use of the wave 
breaking occuring in wake waves generated in a subcrit- 
ical plasma by an ultrashort laser pulse (further called 
"driver"). In a strongly nonlinear wake wave, the elec- 
tron density is modulated in such a way that the electrons 
form relatively thin layers moving with the velocity v p h- 
A counter-propagating laser pulse (called "signal") in- 
teracts with the plasma, perturbed by the driver. Under 
certain conditions, the signal is partially reflected from 
the wake wave, which thereby plays the role of a rcla- 
tivistic mirror. As the wake wave amplitude approaches 
the threshold for wave breaking (i.e., when the electron 
velocity in the wave approaches its wave phase velocity), 
the electric field profile in the wave steepens and, cor- 
respondingly, the singularity is formed in the electron 
density profile. It is important that the singularity af- 
fords partial reflection of the signal with some reflection 
coefficient R. As shown in Appendix B, in general, the 
singularity can be represented as n(X) ex (X — X TO )~ Q , 
1/2 < a < 1. Although the singularity is integrable for 
a < 1, it breaks the geometric optics approximation and 
leads to the efficient reflection of a portion of the radi- 
ation accompanied by the upshifting of the frequency of 
the reflected pulse. When a — » 1, the modulation of the 
electron density is so strong that the reflection coefficient 
becomes to be of the order of unity. 

In order to calculate the reflection coefficient, we con- 
sider the interaction of an electromagnetic wave, repre- 
senting the signal, with a singularity of the electron den- 
sity formed in a breaking Langmuir wave. The electro- 
magnetic wave, given by the z-component of the vector 
potential A z (x,y,t), is described by the wave equation 

d tt A z - c 2 (d xx A z + d yy A z ) + u 2 pe {x - v ph t)A z = 0, (61) 

where lo 2 {X) = 47re 2 n(X)/m e 7 and the electron Lorentz 
factor, 7, is equal to jph at the point where the density, 
n(A), is maximum. 



the form 



with 



d 2 a(0 



\C\ al 



a (0 = z~2 ex P(-^( w '*' - k yV)) 



m e c* 



and 



Here 



n — d h 2 ~ a 'v a ~ 1 



(lu'/cY - k 2 > 0, 



(63) 



(64) 



(65) 



(66) 



and £ = (x — Vphi)7ph, t', k',u)' are the coordinates and 
time and the wave vector and frequency in the boosted 
frame. According to Eqs. (60 1, (651, the coefficient g a 
for a — 2/3 is equal to 



ff2/3 « (2/9)^(1 + ^)l/ 6fc 4/3 7 2/3 i 



(67) 



i.e. G 2/3 = (2/9) V3 (1 + ^)1/6^. 

WKB-approximation gives asymptotic solutions of Eq. 
(63 1 for large \s\ and \(\ 

o(C) « IUi(0{re lsW(c) +pe- JsH/(c) }, (68) 



W(C) 



Co 



(J a 



■s 2 \(\' 



;d£ 



Wi(C) = 1 



2| C | a 



-1/4 



(69) 

(70) 



where r and p are constant. 

We are interested in solutions such that for ( — ► +oo, 
r represents the amplitude of the incident wave (assumed 
to be equal to 1) and p is the reflected wave ampli- 
tude, while in the opposite limit, £ — > — oo, r stands 
for the transmitted wave and p vanishes. Consequently, 
we write: 



A. The case n(X) oc X~ a with a < 1 

Here we consider the general case of the singularity 
formed in the electron density at the point of the wake 
wave breaking, which is defined by 



n(X) 



n G a 
(ft P X)° ' 



(62) 



where G a — const is the dimensionless constant, 1/2 < 
a < 1. 

In the boosted reference frame, moving with the phase 
velocity of the Langmuir wave, the equation (61 1 takes 



r(+oo) = l, |p(+oo)| 5 



R, 



(71) 



|t(-cx^| 2 = T, p(-oo)=0, (72) 

where we assume that the reflection is small, 

R -y and T -» 1 at s -> oo. (73) 



We multiply Eq. (63 1 by da(Q/d£, integrate it over £ 
and obtain 



(^) 2 + »v, 



<2 



<2 



& Ci 



9a da 2 (() 

K\ a d( 



d(. (74) 



According to Eq. (68 1, at large \s\ and |£| we have 



9a 

\c\ a 



T 2 e 2isW(C) 



p 2 e -2tsW(0 



(75) 



Using Eq. ([71]), ([72]), we obtain 

2 



R a — 



1 
Is~ 2 



(«) + S V«> 



g^ T 1 rfa 2 (Q 

4 S 2 y ici« c?c 



d£. 



(76) 



Where the main term in the integrand in the right hand 
side of Eq. d76| for large s is 



(77) 



1 da 2 (Q) 2isW^{C,)ey^{2isW{0) 



IC| Q d( 



2is 



\c\ a 

exp [ 2is£ - 



ig a( 



l-Q 



(l-a)s 



This approximation fails in an interval around £ = 0. 
The interval size decreases at large s and small g a thus 
the intergration over this interval yields negligibly small 
contribution to the reflection coefficient. Neglecting s _1 - 
tcrm in the exponent in Eq. (77), we obtain 



R, 



ig a T(l - a) . /7ra\ 



(78) 



where T is the Gamma function 
coefficient is 



v 2 y 

For a = 1/2, this 



Rl/2 — 



rl/2 



Sl/2 



2s 3 / 2 



For a = 2/3, this coefficient is 

i 3 1 /2r(i/ 3 ) 52/3 



R2/3 — 



(2s) 4 /3 



(79) 



(80) 



For a — ► 1 and fixed s and <7 Q , the main term in the 
expansion of the coefficient (78 1 over (I — a) is R Q « 
ig<x8~ (1 — a)^ 1 . The above used approximation implies 
a smallness of the relection coefficient compared to unity, 
i.e. the condition s S> (1 — a) -1 must be fulfilled. In 
the case a = 1, which requires special consideration, this 
condition can not be satisfied. The case a = 1 is analysed 
below. 



B. The case ra(X) oc X -1 

Here we consider the problem of EMW reflection at 
the following electron density profile 



n(X) 



noGi 
*p|X| 



(81) 



Since this singularity is non-integrable, the thin layer ap- 
proximation is not applicable. As stated above, this pro- 
file formally requires an infinite number of electrons and 
its formation takes an infinite time. We use the standard 



matching technique. Substituting a = 1 to Eq. (63 1 we 
obtain the equation 



dMO 



9i 



o(C) = 



(82) 



d( 2 V ICI. 

where £ is the coordinate in the boosted frame of ref- 



erence and gi is defined by Eqs. (65 The solution to 



this equation can be expressed in terms of the conflu- 
ent hypergeometric functions of the first, M{a, b 1 z), and 
the second kind, U(a, b, z), also known as the Kummcr's 
function of the first and second kind 33 1 : 



-is£ 



o+(C) = it 

+4(M ( 1 
in the region £ > and 



g 1 T(ig 1 /2s) 

igi 



u 



m 

2s 



, 0, 2is( 



2 s 



2,2is( 



(83) 



o-(C) 



-isQ 



3i r (-*5i/2s) V 2s 



-C2CM 1 



m 

2s 



,2,2is( 



(84) 



for the negative coordinate £ < 0. These functions should 
be equal to each other at £ = for all s, i.e. 



o+(+0)-o_(-0) = 



^^MfH* <*> 



which entails c x = c{ . The derivatives of a+(£) and 
a_(t) with respect to £ should differ at £ = for all s by 



+£ 



a',(+0)-a'_(-0) = V.P. 



«(0 

ICI 



d( = 



2ct 



— ln(e) sinh ( — 

7T W V 2 



g± 

2.s 



(86) 



It yields 



=+- C +sinh(^)|1 



id(s)} , 



where 



d(s) 



47E-M 



2 log (2s) 



«5i 
2s 



+ V 1 + 



2s 



(87) 



is real monotonously growing function for g\ > and 
s > 0, 7e-m = 0.577... is the Euler-Mascheroni constant 
and ip(z) is the digamma function |33j . 
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The series expansions of both the branches, (83) and 
(84 1, at £ — > ±oo contain the 
exp (isQ and exp(— isQ. The 



responds to the incident and tr 
the negative exponent - to the 
ously there is no reflected wave 
at £ — > — oo, i.e. the term in 
exp (— is() should vanish. As a 
tionship between cj" and cj~: 



terms proportional to 
positive exponent cor- 
ansmitted wave, while 
reflected wave. Obvi- 
bchind the singularity 
a_(£) proportional to 
result we find a rela- 



cosh 



*(£)+«->(£)<<«] ■ < 88 > 



Now we can find the amplitudes of the reflected, p, and 
transmitted, r, waves as the ratios of the reflected-to- 
incident and transmittcd-to-incident waves, respectively: 




P(s) 



_ T(-i gi /2s) igi/a 



T(iffi/2a) 



(2s) 

exp(-7T5i/2s) 



t(s) 



cosh(ngi/2s) 
T(-i gi /2s) 



T(i gi /2s) 

CXp(-7T5l/s) 



+ isinh(7r<7i/2s)d(s) 
(2s) l9l/s 



cosh(7T5i/2s) + ismh(irgi/2s)d(s) 

The reflection coefficient is 

1 



, (89) 



(90) 



FIG. 5: Coefficients of reflection (Ri = p ), transmission 
(Ti = t 2 ), and absorbtion (Ai = 1 — Ri — Ti) as functions of 
the wave number s for the case n(X) oc X -1 . 



There is a way to increase the number of reflected pho- 
tons by using multiple wake wave periods with the elec- 
tron density singularity in each period. Such the con- 
figuration is typically formed behind the ultrashort laser 
pulse in the underdense plasma. It is easy to show that 
when the pulse interacts with a train of wake waves, the 
number of reflected photons is given by the sum 



Rx = \ P \ 2 = i- 



coshW 2s ) + d( 5 ) 2 sinhV gi / 2s r Nr ' ph = ^B 1 - R r - ».** [i - (i - R)"1 



the transmission coefficient is 

2 exp(-27rgi/s) 



Ti 



cosh^(7r 5l /2s) + d(s) 2 smh 2 (i:g 1 /2s) 



(91) 



(92) 



Here M w is a number of the wake wave periods. Then 
the energy of the reflected electromagnetic radiation is 



f — f 



We note that for the density profile n oc |£| _1 the sum 
of the reflection and transmission coefficients is less than 
one, Ri + Ti = \p\ 2 + \t\ 2 < 1. This corresponds to a 
partial absorbtion of the wave energy in the region of the 
singularity, similarly to the case n ex C _1 anylized in Rcfs. 
[34] and [35]. The dependence of the reflection, transmis- 
sion and absorption, Ai = 1 — Ri — Ti, coefficients on the 
wave number, s, is shown on Fig. [5] 



IV. LASER PULSE INTERACTION WITH 
MULTIPLE RELATIVISTIC MIRRORS 

A number of the photons reflected at the semi- 
transparent mirror depends on the regime of nonlinear 
wake wave breaking. In Refs. [3J El [2J [T^] and above it 
is shown that although the reflection coefficient is small, 
it is not exponentially smal. For example, in the case of 
typical singularity with n(X) oc X~ 2//3 the reflection co- 
efficient in the boosted frame of reference is R s=s I/27 
i.e. a number of photons reflected at a single breaking 



1+0. 

1-/3. 



M 



M 



[1 - (1 - R)"»] 



(94) 



In the limit Khf w < 1, we have N r ^ p h ~ N s ^ p hRAf w , with 
the reflected electromagnetic energy given by 



In the limit Af w 
to 



£r w (4RA^/ 7 m)£ s . (95) 

-¥ 00, the reflected back energy equals 



t r — 4j M c s . 



(96) 



3 



wake is 5N, 



pi 1 



N s ph/2'yf iI . Here N s p h is a number of 



photons in the incident (signal) pulse. 



We see that the energy of the reflected pulse can be many 
times greater than that of the incident pulse. The energy 
gain is due to the momentum transfer from the wake wave 
to the reflected radiation. 

Under conditions which are expected to be realized in 
a moderate intensity laser — gas target interaction, the 
plasma inhomogeneity can play an important role. If 
the driver pulse normalized amplitude, a — eE/m e toc, 
is not large enough to excite the breaking wake wave, 
i.e. a 2 < 2j p h, the wave reaches the breaking threshold 
due to growth of its wave number. As well known, the 
Langmuir waves belong to a class of the waves with con- 
tinuous frequiency spectrum. In inhomogeneous medium 
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the wave vector, k p , and wave frequency, u> p , are related 
to each other via the equation 



<9 t k E 



-Vw„ 



(97) 



The frequency dependence on the coordinates can also be 
due its dependence on the wave amplitude in the case of 
transversally inhomogeneous laser pulse driver, which is 
typical for the laser plasma interaction. As a result, the 
wake wave breaks when w p /|k p | « v e , i.e. when the wave 
phase velocity, u>p/|k p |, becomes equal to the electron 
quiver velocity v e (for details see Refs. [33J ). In this 
regime the breaking wake wave gamma factor, j p h, is 
determined not by the plasma density, but by the driver 
laser pulse amplitude, and it is equal to 7 e . 



V. CONCLUSION 

The plasma wave breaking leads to the formation of 
moving caustics in the plasma flow, which correspond 
to singularities of the particles distribution in the phase 
space. The singularity structure depends of the break- 
ing conditions determined by the parameters of plasma 
and the driver laser pulse which excites the plasma wave. 
We found the structure of typical singularities appearing 
in the plasma particle density during the wave breaking. 
The singularity in the electron density, moving along with 
the wake wave, acts as a flying rclativistic mirror for the 
counter-propagating electromagnetic radiation, reflecting 
the light and up-shifting its frequency. Such the reflection 
results in generation of an ultra-intense electromagnetic 
radiation. We found the coefficients of the reflection of an 
electromagnetic wave at the singularities (crests) of the 
electron density piled up in the most typical regimes of 
strongly nonlinear wave breaking in collisionlcss plasmas. 
Although the reflection coefficients are small, they are 
not exponentially small, as it is usual in the geometric- 
optic approximation of the over-barrier reflection from 
the electron density modulations, applicable for the wake 
wave far below the breaking threshold. The efficiency of 
the photon reflection can be substatially increased by us- 
ing the regimes providing high order sigularity formation 
(see Eq. (78 1 or by realizing the laser pulse reflection 



from several subsequent density singularities, which can 
be naturally produced in the long enough wake wave. 
We show that in the multiple-mirror configuration, the 
energy of the counter-propagating laser pulse reflected 
back and pumped by the wake wave can be increased by 
the factor of &y%[ ■ This paves a way towards producing 
high efficiency compact sources of hard electromagnetic 
radiation. 
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Appendix A 

Plasma flow breaking in the ponderomotive ion 
acceleration 

Here we consider another important example when 
one can neglect the space charge effects. The quasineu- 
tral plasma dynamics is provided by the electromagnetic 
wave interaction with the electron-ion plasmas in the 
limit of extremely high light intensity in the so-called 
"radiation pressure dominated regime" (for details see 
Refs. [T3J mi US]). In the quasineutal plasma motion, 
the longitudinal components of the electron and ion ve- 
locities are equal to each other, i. e. the longitudinal 
component of the electron momentum is much less than 
the ion momentum. Due to conservation of the general- 
ized momentum, the transverse components of the elec- 
tron and ion momenta are determined by the interaction 
with the electromagnetic wave and are equal to m e ca(X) 
and rriiC(Mi(X), respectively, where /i = m e /m i is the 
electron-to-ion mass ratio. It is easy to show that in this 



model instead of the integral ( 26 ) we have the following 
conservation law 



M— + — 1^ ph = M + #- (A.l) 
m e c rriiC; 



Ji + A«7e 



Here p e , 7 e and p t , 7; = \J 1 + (p,a(X)) + (pi/rriiC) are 
momentum and the Lorentz factor for the electron and 



ion, respectively. In the right hand side part of Eq. (A.l I 
the constant H is given by H — 7^0 — (pio/rriic) /3 p h- 

Using the quasineutrality condition, Pi/m^i = 
p e /m e "f e , and smallness of the parameter /1, we reduce 
Eq. (|A.1| to 




TUiC I miC 



/3 ph = 1 + H (l - Vl + a 2 (X)) , 

(A.2) 



and find the ion momentum 



H(X)0 ph -JH*(X)-l 
pi = rmc ; • (A.3) 



pii 



Here H(X) is a function given by the expression 
1 + ix (l - v/l + a 2 (X)) 



H(X) = 



1-A 2 



(A.4) 



P h 
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If th e laser pulse amplitude is so large that a > 
^ _1 -^/(l + /i) 2 — fi 2 , the function H(X) vanishes at 
X = X;, r and the ion momentum and energy become 

equal to TOiC/? ph / ' Jl - /3^ h and rrnc 2 j J\ - (3 2 



J P h- 



The 



plasma layer bounces at this point and then moves in 
the same direction as the electromagnetic pulse. Its mo- 



H>r 



to 



mentum increases from TOjC/Sph/ \ /l — P 2 h at X = X 

mjCy(i + /3 p h)/(l — /3 P h) at X— > oo. More detailed anal- 
ysis of this regime of the ion acceleration will be presented 

elsewhere. 



Appendix B 

Geometric properties of caustics in plasma flows 

We found above three kinds of the singularity formed 
due to the nonlinear wake wave breaking: the particle 
density depends on the coordinate as n(X) ex X~ Q , where 
the exponent a is equal to 1/2, 2/3, and 1. This expo- 
nent represents the measure of the singularity strength. 
Here we present a generic description of nonlinear wave 
breaking, based on the geometric properties of caustics 
in the plasma flow. It is convenient to use the Lagrange 
variables, (x 0} t) , related to the Euler coordinates (x,t) 
as 



X = x +S > (x ,t) 



(B.l) 



with £(xo, t) being the electron displacement and v = d t (, 
. As known, the solution to the continuity equation in 
the Lagrange coordinates is given by the expression 



n(x ,t) = 



»o 



J(x ,t) 



for the particle density, where 



J(x ,t) 



d[x,t) 



d(x ,t) 



\l + d x M 



(B.2) 



(B.3) 



is the Jacobian of transformation from the Lagrange 
to the Euler coordinates. The transformation has a 
singularity at the point where the Jacobian vanishes, 
J{xQ,t) — > , which is equivalent to the condition 
dx £, — > — 1 . This singularity corresponds to the wave 
breaking, when the electron density and the velocity (mo- 
mentum) gradient become infinite. 

In the vicinity of the singularity, where the Lagrange 
coordinate is equal to xq — Xbr , i-e for xq = Xbr + Sxq 
with 5xq <C 1 , we can expand the dependence of the 
Euler coordinate on Xbr and Sxq over Sxq 



Xbr 



X = Xfr 

'HI 

gfer 

2 



Sxn 



Sx n 



£,(Xbr - 



Sx ) 

Sxi - 



(B.4) 



Here the subscript 'br' means that the value of a function 
is taken at xq = Xbr , the prime denotes differentiation 



with respect to xq , and we take into account that 1 + 
£ br = . Expanding the Jacobian J(xo,t) in series of 



Sxq we obtain 

J=\C\Sx - 



!<J) 



(J "1)1 



Sxj, 



(B.5) 



If the second derivative ^' r does not vanish, we obtain 
the following relationship between the Lagrange and Eu- 
ler coordinates: Sx ss (2(x — Xbr)/^) 1 ^ 2 - For the elec- 
tron density we have 



77 o 



(2C) 1/2 



yX ^br) 



(B.6) 



This case corresponds to the regime described by Eq. 
d57l. 



In the case when the second derivative of the dis- 
placement vanishes at the singularity, ££' r = 0, and the 
third derivative is nonzero, £' b " ^ 0, we find Sxq ~ 
(6(.T-^ br )/C) 1/3 and 



JlQ 



2 



1/3 



(X-Xbr)- 2 / 3 . 



(B.7) 



In this case we recover the result obtained above: the 
singularity n ex AX~ ' in the wake wave breaking de- 
scribed by Eqs. ( p|60 |. 

If all the derivatives of the displacement ^ with respect 
to the variable xq of the order less than m vanish, we 
obtain from Eq. (B.4 1 the relationship between x and 
Sx a : 



Xbr 



Jm) 
^>br s„m 



It yields for the Jacobian 



J = m 



'Am)' 



1/r, 



\X Xbr) 



and for the electron density 



1 



n - 



Jm) 
\ ^br 



\X Xbr) 



(B. 



(B-9) 



(B.10) 



where the exponent a m is equal to 

m— 1 ,„ 

a m = B.ll 

m 

with m > 2. 

For m> 1 we use the Stirling formula for the factorial 
to! = ^/27^TO m+1/2 exp(-TO + r m ) with l/(12m + 1) < 
r m < 1/12to. Representing the derivative Q™' as ££™ = 
constA"? 1 , where X w f is the wake field wavelength, we 
find from Eq. (B.10 ) that in the limit to, — > oo the density 



distribution is described by 



n 



-l Tlo\ w f 



(B.12) 



13 



Here e = 2.718... is base of the natural logarithm (the 
Napier constant). This dependence of the electron den- 



sity on the coordinate (a: — Xbr) corresponds to Eq. (33 1 



[1] K. Landecker, Phys. Rev. 86, 852 (1952); L. A. Ostro- 

vskii, Sov. Phys. Usp. 18, 452 (1976); V. L. Granatstein, 

et al. Phys. Rev. A 14, 1194.(1967); J. A. Pasour, V. L. 

Granatstein, and R. K. Parker, Phys. Rev. A 16, 2441 

(1977). 
[2] C. Wilks, et al. Phys. Rev. Lett. 62, 2600 (1989); S. V. 

Bulanov and A. S. Sakharov, JETP Lett. 54, 203(1991); 

C. W. Siders, et al. Phys. Rev. Lett. 76, 3570 (1996); J. 

M. Dias, et al. Phys. Rev. Lett. 78, 4773 (1997); J. T. 

Mendonca, Photon Acceleration in Plasmas (Inst. Phys. 

Publ., Bristol, 2001); C. D. Murphy, et al. Phys. Plasmas 

13, 033108 (2006). 
[3] S. V. Bulanov, et al., Kratk. Soobshch. Fiz. ANSSSR 6, 

9 (1991), in Russian; S. V. Bulanov, et al. in: Reviews 

of Plasma Physics. Volume: 22, edited by V. D. Shafra- 

nov, (Kluwer Academic / Plenum Publishers, New York, 

2001), p. 227. 
[4] S. V. Bulanov, T. Zh. Esirkepov and T. Tajima, Phys. 

Rev. Lett. 91, 085001 (2003). 
[5] M. Kando, et al. Phys. Rev. Lett. 99, 135001 (2007). 
[6] A. S. Pirozhkov, et al. Plasma Phys. 14, 123106 (2007). 
[7] V. I. Semenova, Sov. Radiophys. Quantum Electron. 10, 

599 (1967); M. Lampe, E. Ott and J.H. Walker. Phys. 

Fluids 21, 42(1978); W. B. Mori, Phys. Rev. A 44, 5118 

(1991); R. L. Savage, Jr. et al., Phys. Rev. Lett. 68, 946 

(1992); W. Yu, et al. J. Phys. D 26, 2093 (1993); X. Feng 

and S. Lee, Optics Communications 136, 385 (1997). 
[8] S. S. Bulanov, et al. Phys. Rev. E 73, 036408 (2006). 
[9] S. V. Bulanov, N. M. Naumova and F. Pegoraro, Phys. 

Plasmas 1, 745 (1994); R. Lichters, J. Meyer-ter-Vehn 

and A. M. Pukhov, Phys. Plasmas 3, 3425 (1996); V. 

A. Vshivkov, et al. Phys. Plasmas 5, 2727 (1998); A. S. 

Pirozhkov, et al. Phys. Plasmas 13, 013107 (2006). 
[10] V. V. Kulagin, et al. Phys. Plasmas 14, 113101 (2007). 
[11] N. M. Naumova, et al. Phys. Rev. Lett. 92, 063902 

(2004); N. M. Naumova, J. Nees, and G. Mourou, Phys. 

Plasmas 12, 056707 (2005); N. M. Naumova, et al. New 

J. Phys. 10, 025022 (2008). 
[12] V. Yanovsky, et al. Optics Express 16, 2109 (2008). 
[13] Ya, B. Zel'dovich and A. F. Illarionov, Sov. Phys. JETP 

34, 467 (1971); A. D. Steiger and C. H. Woods, Phys. 

Rev. E 5, 1467 (1972); Ya. B. Zel'dovich, Sov. Phys. Usp. 

18, 79 (1975); A. G. Zhidkov, et al. Phys. Rev. Lett. 

88, 185002 (2002); J. Koga, T. Zh. Esirkepov and S. V. 

Bulanov, Phys. Plasmas 12, 093106 (2005). 
[14] S. V. Bulanov, et al. Plasma Phys. Rep. 30, 221 (2004). 
[15] W. Heisenberg and H. Z. Euler, Z. Phys. 98, 714 (1936); 

J. Schwinger, Phys. Rev. 82, 664 (1951); E. Brezin and 

C. Itzykson, Phys. Rev. D 2, 1191 (1970); N. B. Narozhny 

and A. I. Nikishov, Sov. Phys. JETP 38, 427 (1973); V. 

S. Popov, JETP 94, 1057 (2002); N. B. Narozhny, et al. 

Phys. Lett. A 330, 1 (2004). 
[16] A. Mourou, T. Tajima and S. V. Bulanov, Rev. Mod. 

Phys. 78, 309 (2006); M. Marklund and P. K. Shukla, 

Rev. Mod. Phys. 78, 591 (2006); Y. I. Salamin, et al. 

Phys. Rep. 427, 41 (2006). 
[17] S. V. Bulanov, et al. in: "Proceedings of the Interna- 



tional Workshop on Quark Nuclear Physics", Eds.: J. K. 
Ahn, M. Fujiwara, T. Hayakawa, et al. Pusan National 
Univesity Press, PUSAN, 2006, p. 179. 

[18] S. V. Bulanov and A. S. Sakharov, JETP Lett., 54, 203 
(1991); S. V. Bulanov, F. Pegoraro and A. M. Pukhov, 
Phys. Rev. Lett. 74, 710 (1995); Z.-M. Sheng, et al. Phys. 
Rev. E 62, 7258 (2000); N. H. Matlis, et al. Nature Phys. 
2, 749 (2006); A. Maksimchuk et al. Phys. Plasmas 15, 
056703 (2008). 

[19] A. I. Akhiezer and R. V. Polovin, Sov. Phys. JETP 3, 
696 (1956); O. B. Shiryaev, Phys. Plasmas 15, 012308 
(2008). 

[20] A. Einstein, Ann. Phys. (Leipzig) 17, 891 (1905). 

[21] Ya. B. Zel'dovich and Yu. P. Raizer, Physics of Shock 
Waves and High- Temperature Hydrodynamic Phenom- 
ena (Academic, New York, 1967). 

[22] B. B. Kadomtsev, in: Reviews of Plasma Physics. Vol- 
ume: 22, Ed.: V. D. Shafranov, (Kluwer Academic / 
Plenum Publishers, New York, 2001), p. 1. 

[23] S. V. Bulanov, L. M. Kovrizhnykh and A. S. Sakharov, 
Physics Reports 186, 1 (1990); S. V. Bulanov, et al. Phys. 
Rev. E 58, R5257 (1998); G. Lehmann, E. W. Laedke and 
K. H. Spatschek, Phys. Plasmas 14, 103109 (2007); C. 
G. R. Geddes, et al. Phys. Rev. Lett. 100, 215004 (2008) 
A. V. Brantov, et al. Phys. Plasmas, in press (2008). 

[24] T. Esirkepov, et al. Phys. Rev. Lett. 92, 175003 (2004) 
F. Pegoraro and S. V. Bulanov, Phys. Rev. Lett. 99 
065002 (2007); O. Klimo, et al. Phys. Rev. STAB 11 
031301 (2008); X. Q. Yan, et al. Phys. Rev. Lett. 100 
135003 (2008). 

[25] D. Farina and S. V. Bulanov, Plasma Phys. Control. Fu- 
sion 47, A73 (2005). 

[26] M. N. Rosenbluth and C. S. Liu, Phys. Rev. Lett. 29, 
701 (1972). 

[27] T. Tajima and J. Dawson, Phys. Rev. Lett. 43, 267 
(1979). 

[28] L. M. Gorbunov and V. I. Kirsanov, Sov. Phys. JETP 
66, 290 (1987). 

[29] E. Esarey, et al. IEEE Trans. Plasma Sci. 24, 252 (1996). 

[30] P. Sprangle, et al. Optics Communications 124, 69 (1996) 

[31] P. Chen, et al. Phys. Rev. Lett. 54, 693 (1985). 

[32] T. Esirkepov, et al. Phys. Rev. Lett. 96, 014803 (2006). 

[33] M. Abramowitz and I. A. Stegun, Handbook of Math- 
ematical Functions with Formulas, Graphs, and Mathe- 
matical Tables (Dover, New York, 1964) . 

[34] K. G. Budden, Propagation of Radio Waves (Cambridge 
University Press, Cambridge, 1988). 

[35] T. H. Stix, Waves in Plasmas (American Institute of 
Physics, New York, 1992). 



